Contributions of point defects, chemical disorder, 
and thermal vibrations to electronic properties of Cdi_ x Zn^Te alloys 
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We present a first principles study based on density functional theory of thermodynamic and electronic prop- 
erties of the most important intrinsic defects in the semiconductor alloy Cdi-^Zn^Te with x < 0.13. The alloy 
is represented by a set of supercells with disorder on the Cd/Zn sublattice. Defect formation energies as well as 
electronic and optical transition levels are analyzed as a function of composition. We show that defect formation 
energies increase with Zn content with the exception of the neutral Te vacancy. This behavior is qualitatively 
similar to but quantitatively rather different from the effect of volumetric strain on defect properties in pure 
CdTe. Finally, the relative carrier scattering strengths of point defects, alloy disorder, and phonons are obtained. 
It is demonstrated that for realistic defect concentrations carrier mobilities are limited by phonon scattering for 
temperature above approximately 150 K. 

PACS numbers: 61.72.J-, 72.80.Ng, 72.10.Fk, 71.55.Gs 



03 



o 
o 



> 
in 
o 
m 
m 

O 

(N 



X 



I. INTRODUCTION 

Cadmium telluride (CdTe) has found applications in a wide 
range of areas including but not limited to solar cells*! radi- 
ation detectors^ and electro-optical modulators^ CdTe is a 
prime candidate material for room-temperature radiation de- 
tectors thanks to a number of reasons: (1) Its band gap lies 
within the range that maximizes the signal-to-noise ratio by 
balancing the numbers of photogenerated carriers vs thermally 
generated carriers. (2) High purity CdTe crystals exhibit a 
long carrier drift length, characterized by the product of car- 
rier mobility /i and carrier trapping lifetime r4r2 A large /it 
value ensures efficient charge collection and a linear response 
to the energy of incident radiation. (3) Finally, these samples 
possess also high resistivity p, which is required to keep leak- 
age currents at a minimum. 

Further improvements can be accomplished by addition of 
Zn to obtain Cdi-^Zn^Te (CZT) alloys. The band gap grad- 
ually increases with Zn concentration from 1.5 to approxi- 
mately 1.6 eV at x w 0.1, which leads to a further reduction of 
thermal noise and raises the resistivity up to 8 x 10 10 VI cm£2& 

Zn alloying has empirically been found to have additional 
benefits such as oxygen gettering, increased hardness, and re- 
duced degradation of detector performance over time from 
build-up of internal fields under high radiation flux ("polar- 
ization" effect)4- These features have allowed the fabrication 
of CZT gamma detectors with excellent energy resolution bet- 
ter than 1 % at 662 keVA^ The large scale deployment of CZT 
detectors is, however, hampered by high material costs. This 
results from the difficulty to grow large uniform single crys- 
tals, which requires cutting the ingot into individual crystal- 
lites, followed by testing and selection ("crystal harvesting"). 
A variety of defects — from point defects to dislocations and 
secondary phases — limit the performance of the individual 
crystals. Understanding and learning to control the different 
contributions is therefore key to improving the manufacturing 
yield and driving down cost. 

Defects affect resistivity, charge carrier mobilities, and car- 
rier lifetimes in different ways. The resistivity for example is 



largely determined by the concentrations of charged defects, 
which are coupled by the charge neutrality condition and can 
be obtained from the defect formation energiesJ-^ Charge 
carrier mobilities in turn are limited by scattering at lattice 
perturbations and are therefore sensitive not only to point de- 
fects and dislocations but also thermal vibrations. Finally, 
lifetimes are shortened by carrier trapping, which occurs pref- 
erentially at defects that induce deep levels in the band gap 
leading to carrier recombination via the Shockley-Read-Hall 
mechanism. 13 In CdTe and CZT, generally lifetime-limiting 
defects remain the biggest challenge to maximizing device 
performance. 

These effects are modified in the case of alloys. The chem- 
ical disorder introduces a distribution of defect formation en- 
ergies and electronic transition levels, and provides additional 
scattering channels. In this work we compute the formation 
energies and transition levels of native defects in CZT and 
study their dependence on alloy composition as well the pos- 
sible appearance of new defect levels in the gap. Further- 
more, we determine the relative importance of point defects, 
phonons, and alloying for carrier mobilities. This assessment 
is important for developing material optimization strategies. 
For example intentional doping and annealing can be used to 
modify the defect populations with little impact on phonon 
properties. 

There are approaches at different levels of sophistication for 
modeling the thermodynamic properties of an alloy includ- 
ing for example the virtual crystal approximation (VCA)r^i^ 
the coherent potential approximation, 1 - special quasirandom 
structures^ 7 - and alloy cluster expansions^ Here, in order to 
obtain a simple unified approach that allows for the study of 
both the distribution of defect-related thermodynamic quanti- 
ties and scattering rates, we model the alloy disorder by using 
a set of supercells with Zn atoms randomly distributed over 
Cd sites. This approach has been applied successfully by Car- 
valho et cd. to study formation energies of cation-site intrinsic 
defects in CZT— In the following we employ density func- 
tional theory to demonstrate from first principles that (;) point 
defect formation energies increase or stay approximately con- 
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stant as a function of Zn concentration, (it) optical transition 
levels are virtually unchanged upon alloying, and (Hi) mobil- 
ities are limited by phonon scattering at ambient conditions. 
The study focuses on Zn concentrations up to 15%, which is 
the most relevant range for applications in radiation detection. 



H. METHODOLOGY 

A. Point defect thermodynamics 

The temperature dependence of the equilibrium defect con- 
centration obeys the following relation 2 ^ 



cq exp 



AG/ 
k B T 



(1) 



where cq denotes the concentration of possible defect sites and 
AG / is the Gibbs free energy of defect formation, which we 
approximate by the formation energy AEf, assuming that vi- 
brational entropy and pressure-volume term are small^ The 
formation energy for a defect in charge state q is given b y 22,23 

AEf = E de f - £ideal - ^ An,ft + q (E VBM + /i e ) (2) 



where E^f and -Eideai are the energies of the system with and 
without the defect, respectively, An; denotes the number of 
atoms of type i and chemical potential [i\ added to the sys- 
tem, and fi e is the electron chemical potential measured with 
respect to the valence band maximum (VBM), £Vbm- 

Neglecting entropic contributions , the heat of formation for 
the CZT alloy is given by 



AiJ / (CZT, x) = (l- x)fl C d + X^Zn + MTe 

— (\ — T\n hulk — <r;/ bulk — /; buIk 
I 1 x )^Cd x t 1 Zn MTe ■ 



(3) 



Defining /ij = /i bulk + A/i i; where /i bulk is the chemical poten- 
tial of component i in its reference state, we obtain the relation 

Aff/CCZT, x) = (1 - x)A^cd + xA^ Zn + A/ixe- (4) 

Since A/i; has to be non-positive to avoid forming bulk mate- 
rial i, the range of the chemical potentials is constrained by 



< Am < AH f (CZT, x) 



(5) 



We then define the cation-rich limit as (1 — x)A/j,cd + 
xAfizn = and the Te-rich limit by A/ixe = (see Ref. 
l24l) . Given the formation energies at fj, e = for two dif- 
ferent charge states, qi and q2, we compute the equilibrium 
electronic transition level as 

F el (n , _ AE f (q 2 ) - AE f ( qi ) 

91 - q-2 



B. Theoretical carrier mobilities 

In a quasi-classical picture, the motion of free carriers sub- 
jected to a weak electric field is governed by the Boltzmann 



transport equation^ In particular, the electrical conductiv- 
ity and carrier mobilities are related to scattering processes 
at the microscopic scale.— In an alloy system, three major 
carrier scattering channels are active, namely those involv- 
ing phonons, the inherent randomness of the alloy itself, and 
defects . 14 i 15 i 26 ~ — The relative importance of these scattering 
channels can be assessed from first principles using Fermi's 
Golden rule 
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(flAV^Sief-e.,). 



(7) 



Here i and / denote initial and final states, respectively, and 
AV is the relevant scattering potential. To compare the rel- 
ative importance of each scattering channel, each contribu- 
tion should be treated at the same level of approximation. Al- 
though highly accurate first-principles approaches have been 
developed to tackle each channel separately ^22. the compu- 
tational cost associated with the combinatorial explosion of 
point defects in an alloy currently hinders further exploration. 
Therefore, more cost-effective approaches are desirable. In 
the following we present an approximate, but comprehensive, 
approach that forms the basis for our analysis. 

a. Defect scattering For point defect scattering, we be- 
gin with Fermi's golden rule as shown in (0. The orthogonal- 
ity of the initial and final states ensures that only gradients of 
AV contribute to the integral. Thus, we have earlier defined 
a simpler, but cruder, measure termed the relative scattering 
strength by 2 ^ 



M 



dr |V r (AV0| 



(8) 



Note that this quantity is independent of initial and final state 
wave functions. Despite its simple form, we have shown 
that M 2 from Eq. ([8]) is approximately proportional to the 
more computationally expensive Brillouin-zone average over 
all scattering rates given by Eq. (0.— 

b. Alloy scattering To quantify scattering due to alloy 
disorder, we follow a similar procedure as for defect scatter- 
ing. For a given Zn concentration x we calculate the relative 
scattering strength using the scattering potential AV given by 



AV 



alloy 



Kdloy - (1 - a;)VcdTe - xV 2 



ZnTe i 



(9) 



where all terms are computed for the same supercell size and 
lattice constant. We then average over an ensemble of random 
alloy configurations for each x. Substitution of (1 — x) VcdTe + 
xVznTtt by the potential resulting from a VCA calculation at 
the same concentration did not change the results. 

c. Phonon scattering Here we only consider phonon 
scattering in pure CdTe and assume the results to transfer to 
the CZT alloy. This is justified since the vibrational frequen- 
cies and phonon polarization vectors do not change much for 
small concentrations of Zn in CdTe. Again following the per- 
turbative approach as above, we compute the relative scat- 
tering strength of a thermal population of phonon modes by 
taking the difference between a thermally excited atomic con- 
figuration and the ideal zero temperature configuration. We 
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perform the ensemble average over phonon modes by aver- 
aging over different atomic configurations representing ther- 
mally populated phonons. The configuration space was sam- 
pled in two different ways to separate harmonic and anhar- 
monic contributions. In the first approach, cells were gen- 
erated by superimposing harmonic phonons modes with ran- 
dom phases and amplitudes chosen to yield the correct mean 
square displacement for the temperature of interest. In the sec- 
ond approach, anharmonicity was explicitly accounted for by 
sampling decorrelated configurations from ab initio molecular 
dynamics trajectories. 
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C. Computational details 

Density functional theory (DFT) calculations were car- 
ried out in the local density approximation (LDA) using the 
Vienna ab-initio simulation package (VASP)2£ and ultrasoft 
pseudopotentials^ We employed supercells of 216 atoms 
sampled in reciprocal space with 2x2x2 /c -point grids gen- 
erated using the Monkhorst-Pack scheme. 32 The plane wave 
energy cutoff was set to 277 eV and Gaussian smearing with 
a width of 0. 1 eV was used to determine the occupation num- 
bers. For a given Zn concentration we created 20 replicas of 
the original ideal or defective cell and randomly substituted 
Zn for Cd. The cells were subsequently relaxed until ionic 
forces were less than 30 me V/A. The lattice parameters were 
taken from VCA calculations.— A homogeneous background 
charge was added for calculations including a charged defect 
to ensure charge neutrality of the entire cell. The monopole- 
monopole correction of Makov and Payne was applied to cor- 
rect for the spurious interaction of charged defects^ No band 
gap or potential alignment corrections were performed. 

We note that equilibrium transition levels Ef r are distinct 
from optical transition levels as the former determine 
the equilibrium electron chemical potential at which a change 
in the charge state occurs whereas the latter ones correspond 
to optical ("vertical") transitions without the involvement of 
lattice relaxations. Neglecting excitonic or multiplet effects 
we can in principle calculate optical transition levels from to- 
tal energy differences of a defect in charge state q and q ± 1 
in the same geometry. Thus, we need an additional separate 
calculation of the total energy for each relaxed defect geom- 
etry. Due to the very large number of defect configurations 
we here instead make the additional approximation that 
may be calculated by Kohn-Sham eigenvalue differences for 
charge state q. We note that this corresponds to the neglect of 
changes in the DFT double counting term (see Ref. l2lh . 

In this work, we focus on the most important defects as 
determined by earlier studies^ 3 ^ These are the Cd and Te va- 
cancies, Te on Cd antisite (in local Td and C^ v symmetries), 
Te-tetrahedrally coordinated Cd interstitials, and the most im- 
portant Te-Te and Cd-Te dumbbell split interstitials. We refer 
to Ref. 36 for details of the respective geometries. 



FIG. 1 . (Color online) Free energy of mixing of CZT as a function 
of Zn concentration and temperature. Symbols denote the enthalpies 
of mixing calculated using the total energies from the supercell ap- 
proach. The red solid line is the result of a quadratic polynomial fit to 
the lowest energy configurations within the gray area (< 13%). The 
remaining lines represent the free energy of mixing up to 400 K, as 
generated by adding the entropy contribution for a random solution. 



III. RESULTS AND DISCUSSIONS 

A. Ideal CZT alloy 

The adequacy of the supercell approach, outlined in Sect. 
| to mimic the behavior of the ideal alloy was assessed 
by a study of its thermodynamic properties below the melt- 
ing temperature. Firstly, the zero Kelvin enthalpy of mixing 
was obtained by fitting a second order polynomial to the low- 
est energy configurations for Zn concentrations up to 13%, 
as displayed in Fig. Q] This corresponds to the enthalpy of 
mixing for pairwise interaction energies and a random dis- 
tribution of atoms on the sub-lattices. 37 Thus, the system is 
completely immiscible at absolute zero. The critical temper- 
ature (T c ) was then approximated by assuming a completely 
random solution, yielding T c = 340 K. In contrast, the VCA 
(results not shown) erroneously predicts a miscible solution at 
all temperatures. 

The literature is unfortunately riddled with disparate esti- 
mates of T c . Model calculations based on the CALPHAD 
approach 38 give 701 K, while cluster expansion calculations 
based on first principles energies^ predict 605 K. Experi- 
mental estimates also vary, with one careful study report- 
ing T c = 435 K from a fit to extended x-ray absorption 
fine-structure characterization of phase separating alloys, but 
also T c = 340 K by extracting parameters from the liquid- 
solid phase diagram4S Earlier experimental studies had re- 
ported and accepted higher values for T c , but it is pointed out 
that slow diffusion at low temperatures makes the determina- 
tion difficult by observation of phase separation, which might 
grossly overestimate the binodal stability temperature. 

In addition, we calculated the variation in band gap as a 
function of Zn concentration and compared to experimental 
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FIG. 2. (Color online) Deviation from linear dependence of band 
gap as a function of Zn concentration. The solid yellow line and the 
red open circles denote the VCA and supercell results, respectively. 
The calculated values are shifted with a linear function to coincide 
with the room-temperature experimental band gap of CdTe and ZnTe 
at x = and x = 1, respectively. Experimental values recorded at 
room-temperature are shown with blue symbols, see Refs. Rll - B^ . 
Data recorded below 10 K arc shown with green symbols, sec Refs. 
|46| andl47l. The latter have been shifted by a constant to coincide with 
the room-temperature data. 



data, as shown in Fig. [2] The randomized supercell approach 
accurately reproduces the experimental trend, while the VCA 
exhibits incorrect bowing. In fact, even the sign of the bowing 
is incorrect with VCA over almost the entire range. 



B. Formation energies and electronic transition levels 

We display in Fig.[3]the formation energies (lines) and elec- 
tronic transition energies (symbols) as a function of the elec- 
tron chemical potential under Te-rich conditions for the spe- 
cific case of x = 0.08. To summarize, the most dominant 
defects are the double acceptor V^? d ~ , double donor Cd^j e , 

and the electrically neutral Te^ d . In fact, the relative order- 
ing in terms of formation energies and electronic transition 
levels for the different defects is unchanged with respect to 
CdTe.— However, in the case of CZT a few defects, most no- 
tably Vje and Cd^xe, display a strong variation of formation 
energies and electronic transition levels between different al- 
loy configurations (supercells). This effect is highlighted in 
Fig. |4j where AEf and Ef r are plotted separately for each 
configuration as a function of Zn concentration. 

In principle, strong variations in the formation energies are 
a tell-tale sign of defects being sensitive to local structure vari- 
ations. The impact of Zn composition and alloy disorder on 
formation energies will now be discussed the cases of Vj e and 
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FIG. 3. (Color online) Dependence of the formation energies on the 
electron chemical potential under Te-rich conditions for the most im- 
portant intrinsic defects in Cdi-^Zn^Te (x — 0.08). Each line corre- 
sponds to a separate random Zn configuration. The electronic (charge 
state) transition energies are denoted by the symbols. 



1. Tellurium vacancy 



In contrast to the anion vacancies in III-V 
semiconductors,^2£^Il the Vr e defect is not subject to a 
Jahn-Teller distortion.- 1 Rather, in an unrelaxed supercell 
containing a single neutral Te vacancy, the a\ defect level 
resides inside the band gap and is fully occupied, whereas 
the triply degenerate ti level is completely empty. Thus no 
energy can be gained by breaking the Td symmetry. When 
allowing the ions to relax, the neighboring Cd ions move 
symmetrically inward by a distance of 2.0 A toward the 
vacant site for q = and outward by 3.5 A for q = 2. The 
inward relaxation for q = causes the a\ defect level to 
lower its energy and hybridize with valence band states. In 
CZT the Vt? + geometry is virtually unchanged with respect 
to the defective CZT, whereas the geometry of the q = state 
displays drastic unsymmetric variations. In particular, the 
nearest neighbor (NN) Zn atoms tend to not relax inwards 
as the Cd NNs do. We however emphasize that in some 
configurations Cd NN also display this behavior. As shown 
Fig.[5k, this manifests itself as two two distinct average Vj e - 
NN distances (d< and d>). Furthermore, the unsymmetric 
configurations (d>), which in most cases involves a Zn NN, 
tend to have lower formation energies. Note that this effect is 
not observed for charge state q = 2. However, as the overlap 
of the formation energies for the two distinct groupings is 
appreciable, the large variation cannot be correlated to the 
number of NN Zn neighbors or local geometry alone. 

Figure displays the variation of formation energies with 
respect mean square displacement J2i °f a ^ atoms. It can 
be seen that the formation energies decrease with increasing 
amount of relaxation. The reduction is roughly linear for the 
geometries that maintain the local Td symmetry. Thus we con- 
clude that the variation in formation energies is due to atomic 
relaxation outside the nearest neighbor shell. This is consis- 
tent with the earlier observation that the occupied a\ defect 
state hybridizes with extended valence band states. 
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FIG. 4. (Color online) (a) For each defect and charge state, the formation energy at the VBM for Te-rich conditions is plotted for the set of 
configurations for five different concentrations ranging from 5.5 to 13.0% (left to right, indicated by arrow), (b) For each defect, the electronic 
transition levels are plotted for the set of configurations for five different concentrations ranging from 5.5 to 13.0% (left to right, indicated by 
arrow). 
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FIG. 5. (Color online) (a) Vt c formation energy at the VBM as a 
function of average Vrc-nearest neighbor distance. The Zn concen- 
tration is indicated by the different symbols. Red and green sym- 
bols correspond to charge states q = and q — 2, respectively. 
Dashed line indicate separation of configurations with small (d<) 
and large(d>) Vtc-NN distances, (b) Vtc formation energy for q = 
at the VBM as a function of displacements of all atoms with respect 
to the defect free system. Again, symbols denote Zn concentration 
and colors indicate small (d<) and large(d>) Vtc-NN distances. The 
line is drawn to guide the eye. 



2. Cadmium interstitial 



FIG. 6. (Color online) Formation energy at the VBM as a function 
of average Cdijc -nearest neighbor distance. The symbols and colors 
indicate the Zn concentration and the number of Zn atoms in the 
next-nearest neighbor shell. 



of average NN distance and number of next-NN Zn atoms. 
One finds the average formation energy to generally increase 
with Zn concentration (see Fig. @]). Furthermore, the large 
variation of formation energies is found to be strongly depen- 
dent on the number of next-NN Zn atoms which also corre- 
lates with the average NN distance. 



The analysis of the Cd; Te geometries is more straightfor- 
ward. In CdTe, the Cd?^ e is tetrahedrally coordinated to four 
Te, each at a distance of 2.85 A. By random substitution of 
Cd by Zn, the formation energies spread to a range of about 
0.35 eV. In Fig.|6]we show the formation energies as a function 



3. Alloying vs strain 

The general trends for all defects studied here are clear 
across the concentration range of interest: the formation en- 
ergies increase or stay constant as a function of Zn concen- 
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tration for all defects and charge states, except Vcd- Carvalho 
et al. examined the correlation between point defect forma- 
tion energies for cation-site defects in CZT and volumetrically 
strained CdTe and concluded that in the range of a: < 0.5, the 
volume change induced by alloying is the dominant factor in 
the change of formation energies.—' Thus, they posit that accu- 
rate defect formation energies for the alloy may be predicted 
by performing calculations with the pure binary compound 
using the lattice constant of the alloy. To the benefit of this 
conclusion, we show in Fig.|7]that usually point defect forma- 
tion energies in CZT and strained CdTe indeed exhibit some 
general overall correlation. However, the absolute values of 
formation energies as well as their variation with volume are 
quantitatively rather different. In fact, the only defect showing 
reasonably good quantitative agreement is the Cd interstitial. 
Hence we conclude that formation energies calculated using 
the strained binary compound in general have small, if any, 
predictive power. 



C. Optical transition levels 

The optical transition levels for defects with levels inside or 
near the band gap have been extracted from the Kohn-Sham 
band structure of the defect cell and are displayed in Fig. [8] 
as a function of Zn concentration. The black solid lines indi- 
cate optical transitions in pure CdTe, and each colored line 
represents a defect state in a particular randomized super- 
cell. In general, the average for each group of levels coin- 
cides roughly with the CdTe result. The alloy configurations 
show various spreads for each defect, arising for the same rea- 
sons as the variations in the formation energies described in 
Sect. IHI Bl Furthermore, with a few exceptions (notably, Te^j 
and (Te— Te)j e <110 >X the position and spread of the transi- 
tions are fairly insensitive to Zn concentration. Since no new 
levels are introduced inside the gap (nor do any existing levels 
move appreciably deeper into the gap) from Zn alloying, we 
conclude that the effect of defect-mediated electron-hole re- 
combinations via the Shockley-Read-Hall mechanismi 3 - does 
not increase upon Zn alloying. 



D. Defect, alloy, and phonon scattering 

The relative carrier scattering strengths for the intrinsic de- 
fects as a function of Zn concentration are shown in Fig. [9] 
These scattering strengths correspond to one defect in a 216- 
atom cell and hence correspond to a defect concentration on 
the order of 10 20 cm -3 . Compared to pure CdTe, the scatter- 
ing strengths are virtually unchanged and exhibit little scatter 
with alloy configuration. A large lattice distortion with re- 
spect to the ideal cell gives rise to scattering potentials with 
large gradients and hence strong scattering. In accordance, 
the split interstitials and Te vacancies are by far the strongest 
carrier scattering centers. It is also quite natural to expect that 
scattering strengths of defects in pure CdTe would serve as 
a lower bound since the alloy disorder could introduce more 
complex relaxed geometries as well as stronger scattering po- 



tentials. Indeed, this is in fact observed for all defects in the 
present study with the exception of the neutral Te vacancy. We 
recall, however, that the nearest neighbors of Vj e are subject 
to a large inward relaxation, unless it is a Zn atom (see Sect. 
HUB I ). Hence, the defect structure for the alloy is actually less 
distorted with respect to the ideal case when a Zn atom sits 
next to the vacancy, explaining this apparent anomaly. 

The last column in Fig.[9]compares the relative scattering of 
randomized 216-atom cells with respect to the ideal structure 
at the same lattice constant, representing the effect of alloying. 
In this case the scattering strength is strongly dependent on the 
Zn concentration and increases mono tonic ally. As expected, 
the effect of alloying on carrier scattering far outweighs defect 
scattering for any realistic defect concentration. 

The contribution to carrier scattering from phonons is 
shown in Fig. [TO) comparing the results from MD simula- 
tions to the harmonic approximation. The latter exhibits a 
linear temperature dependence whereas the MD data curves 
upward due to anharmonic effects. The deviation between 
harmonic and MD data points is a measure for the degree of 
anharmonicity, which already at room temperature amounts 
to 20% of the relative scattering strength. Compared to defect 
and alloy scattering, these results predict that electronic trans- 
port is limited by phonon scattering for temperatures above 
approximately 150 K. 



IV. CONCLUSIONS 

The effects of alloying on pure bulk properties, thermody- 
namic properties of intrinsic point defects, and carrier scatter- 
ing rates were examined for CZT alloys with up to 13% Zn. 
The defect formation energies of the most important intrinsic 
defects increase or remain constant with increasing Zn con- 
tent with the exception of the Cd vacancy and to some extent 
the neutral Te vacancy. The latter defect also deviates from 
the standard inward relaxation of nearest neighbor atoms in 
the presence of Zn NNs, creating a very complex and rich 
energy landscape as a function of Zn coordination (reminis- 
cent of the effects of In coordination in the dilute nitride al- 
loy GalnNA a 52 - 53 ). A previously reported relation between 
strained CdTe and CZT— was revisited, revealing that the de- 
pendence of defect formation energies on Zn concentration 
is at most qualitatively similar to the behavior deduced from 
volumetrically strained CdTe. 

Optical transition levels exhibit small variations with Zn 
content and introduce no additional defect levels within the 
gap. Thus Zn alloying should not increase the defect-mediated 
electron-hole recombinations. 

The relative carrier scattering rate concept introduced in 
Ref. |29| was applied here to alloy and phonon scattering as 
well as defect scattering. The intrinsic defect scattering in 
CZT was shown to be very similar to pure CdTe. The scat- 
tering rates in CdTe were shown to serve as as lower bound 
for all defects except the Te vacancy. The latter behavior was 
explained in terms of local lattice relaxation associated with 
the distribution of Zn NNs. We furthermore established the 
relative importance of the different scattering channels and 




FIG. 7. (Color online) Formation energies at the VBM under Te-rich conditions for the 6 most important intrinsic defects in CdTe as a function 
of compressive strain (lines) and CZT as a function of volume (symbols). 
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FIG. 8. (Color online) Kohn-Sham defect levels for the most im- 
portant intrinsic defects in CZT. For each defect, levels of all the 
randomized cells are shown for five different concentrations ranging 
from 5.5 to 13.0 % (left to right, indicated by arrow). Red, blue, and 
green line colors denotes an empty, half-filled, and fully occupied 
level, respectively. The solid black lines show the level positions in 
the pure CdTe case. 



FIG. 9. (Color online) Relative scattering strengths (Eq. [8]l for the 
most important intrinsic defects in CZT. For each defect, scattering 
strengths are shown for five different concentrations ranging from 
5.5 to 13.0 % (left to right, indicated by arrow). The solid black lines 
show the relative scattering strengths for the pure CdTe case. The last 
column (Alloy) shows the effect of alloy scattering at the same level 
of approximation. The data are arbitrarily normalized to the value of 
the strongest intrinsic defect scattering channel. 
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FIG. 10. (Color online) Relative thermal scattering strengths for 
CdTe as a function of temperature, calculated using both a harmonic 
approximation (solid line) and molecular dynamics (dashed line). 



1 J. Britt and C. Ferekides, Appl. Phys. Lett. 62, 2851 (1993). 

2 G. Knoll, Radiation Detection and Measurement, 4th ed. (John 
Wiley & Sons, 2010). 

3 R. Triboulet and R Siffert, CdTe and Related Compounds; 
Physics, Defects, Hetero- and Nano-structures, Crystal Growth, 
Surfaces and Applications: Crystal Growth, Surfaces and Ap- 
plications, European Materials Research Society Series (Elsevier 
Science, 2009). 

4 S. Del Sordo, L. Abbene, E. Caroli, A. M. Mancini, A. Zappettini, 
and P. Ube rtini, Sensors 9, 3491 (2009). 

5 A. Owens, |j~Syhc. Rad. 13, 143 (2006)| 

6 R. Redus, J. Pantazis, T. Pantazis, A. Huber, and B. Cross, 
| IEEE Trans. Nucl. Sci 56, 2524 (2009)1 

7 T. Schlesinger, J. Toney, H. Yoon, E. Lee, B. Brunett, L. Franks, 
and R. James, Mate r. Sci. Eng. R 32, 103 (2001)| 

8 A. Melnikov, A. Sigov, K. Vorotilov, A. Davydov, L. Topalova, 
and N. Zhavoronkov, J. Cryst. Growth 197, 666 (1999). 

9 F. Zhang, Z. He, D. Xu, G. Knoll, D. Wehe, and J. Berry, IEEE 
Trans. Nucl. Sci 51, 2427 (2004). 

10 S. A. Awadalla, H. Chen, J. Mackenzie, P. Lu, K. Iniewski, 
P. Marthandam, R. Redden, G. Bindley, Z. He, F. Zhang, 
M. Groza, A. Burger, D. R. Mayo, and C. L. Sullivan, in Nuclear 
Science Symposium Conference Record, 2008. NSS '08. IEEE 
(2008) pp. 58 -62. 

1 1 P. Erhart and K. Albe, | Journal of Applied Physics 104, 044315 (2008) 

12 P. Erhart, D. A. berg, B. W. Sturm, K.-J. Wu, and V. Lordi, Ap- 
plied Physics Letters 97, 142104 (2010). 

13 W. Shockley and W. T. Read, Phys. Rev. 87, 835 (1952). 

14 L. Nordheim, Ann. Phys. (Berlin) 9, 607 (1931). 

15 L. Nordheim, Ann. Phys. (Berlin) 9, 641 (1931). 

16 P. Soven, Phys. Rev. 156, 809 (1967). 

17 A. Zunger, S.-H. Wei, L. G. Ferreira, and J. E. Bernard, Phys. Rev. 
Lett. 65, 353 (1990); S.-H. Wei, L. G. Ferreira, J. E. Bernard, and 
A. Zunger, Phys. Rev. B 42, 9622 (1990). 

18 D. B. Laks, L. G. Ferreira, S. Froyen, and A. Zunger, Phys. Rev. 
B 46, 12587 (1992). 

19 A. Carvalho, A. K. Tagantsev, S. Oberg, P. R. Briddon, and 
N. Setter, Phys. Rev. B 81, 075215 (2010). 

20 A. R. Allnatt and A. B. Lidiard, Atomic Transport in Solids (Cam- 



bridge University Press, Cambridge, 2003). 

D. Aberg, P. Erhart, A. J. Williamson, and V. Lordi, Phys. Rev. B 

77, 165206 (2008). 

G.-X. Qian, R. M. Martin, and D. J. Chadi, Phys. Rev. B 38, 7649 

(1988). 

S. B. Zhang and J. E. Northrup, |Phys. Rev. Lett. 67, 2339 (1991)1 

L. E. Ramos, J. Furthmiiller, F. Bechstedt, L. M. R. Scolfaro, and 

J. R. Leite, J. Phys.: Condens. Matter 14, 2577 (2002). 

G. Mahan, Many-particle physics, Physics of solids and liquids 

(Kluwer Academic/Plenum Publishers, 2000). 

P. Yu and M. Cardona, Fundamentals of Semiconductors: Physics 

and Materials Properties, Graduate Texts in Physics (Springer, 

2010). 

F. Murphy- Armando and S. Fahy, 
|Phys. Rev. Lett. 97, 096606 (2006)| 

S. Joyce, F. Murphy- Armando, and S. Fahy, 

|Phys. Rev. B 75, 155201 (2007)] 

V. Lordi, P. Erhart, and D. Aberg, 

|Phys. Rev. B 81, 235204 (20T0)j 

G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993); 49, 14251 
(1994); G. Kresse and J. Furthmiiller, 54, 1 1 169 (1996); Comput. 
Mater. Sci. 6, 15 (1996). 

D. Vanderbilt, [Phys. Rev. B 41, 7892 (1 990) 

H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976). 
The VCA for total energies, forces, and stress was implemented 
into VASP using the approach of Bellaiche and Vanderbilt.— Note 
that this implementation is separate from the current official re- 
lease. 

G. Makov and M. C. Payne Phys. Rev. B 51, 4014 (1995). 
S.-H. Wei and S. B. Zhang, |PhysT Rev. B 66, 155211 (2002)| 
P. Erhart, D. Aberg, and V. Lordi, In preparation. 
A.-B. Chen and A. Sher, Semiconductor alloys: physics and ma- 
terials engineering, Microdevices Series (Plenum Press, 1995). 
A. Marbeuf, R. Druilhe, R. Triboulet, and G. Patriarche, 
|J. of Crys. Growth 117, 10 (1992)| 

S.-H. Wei, L. G. Ferreira, and A. Zunger, 

|Phys. Rev. B 41, 8240 (1990)] 

N. Motta, A. Balzarotti, P. Letardi, A. Kisiel, 
M. Czyzyk, M. Zimnal-Starnawska, and M. Podogorny, 



9 



|J. Cryst. Gr owth 72, 205 (1985) 

O. Zelaya-Angel, J. G. Mendoza-Alvarez, M. Becerril, 
H. Navarro-Contreras, and L. Tirado-Mejfa, J. Appl. Phys. 
95, 6284 (2004). 

J. Kimmerle, R. Menner, H. Scock, and A. Valera, Thin Solid 
Films 126, 23 (1985). 

K. Prasada Rao, O. M. Hussain, K. T. R. Reddy, P. S. 
Reddy, S. Uthanna, B . S. Naidu, and P. J. Reddy, 
|Opt. Mater. 5, 63 (1996)| 

V. A. Tyagai, S. O. V., V. N. Bondarenko, N. I. Vitrikhovskii, V. B. 
Popov, and N. Krasiko, Sov. Phys. Solid State 16, 885 (1974). 
K. K. Chattopadhyay, A. Sarkar, S. Chaudhuri, and A. K. Pal, 
Vacuum 42, 1113 (1991). 



F. G. Sanchez-Almazan, H. Navarro-Contreras, G. Ramfrez- 
Flores, M. A. Vidal, O. Zelaya-Angel, M. E. Rodriguez, and 
R. Baquero, J. Appl. Phys. 79, 7713 (1996). 

47 J. L. Reno and E. D. Jones, |Phys. Rev. B 45, 1440 (1992)| 

48 J. Bernholc and S. T. Pantelides, P hys. Rev. B 18, 1780 (1978)] 

49 D. N. Talwar and C. S. Ting, |Phys. Rev. B 25, 2660 (1982)| 

50 M. J. Puska, J. Phys.: Condens. Matter 1, 7347 (1989). 

51 S. Lany, V. Ostheimer, H. Wolf, and T. Wichert, Phys. B: Cond. 
Matt. 308-310, 958 (2001). 

52 V. Lordi, V. Gambin, S. Friedrich, T. Funk, T. Takizawa, K. Uno, 
and J. S. Harris, Phys. Rev. Lett. 90, 145505 (2003). 

53 V. Lordi, H. B. Yuen, S. R. Bank, M. A. Wistey, J. S. Harris, and 
S. Friedrich, Phys. Rev. B 71, 1 25309 (2005). 

54 L. Bellaiche and D. Vanderbilt, [Phys. Rev. B 61, 7877 (2 000) 



